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Abstract. We revisi t experimental studies performed by Ek- 
man on dead-water dEkmani 1 19041) using modem techniques 
in order to present new insights on this peculiar phenomenon. 
We extend its description to more general situations such 
as a three-layer fluid or a linearly stratified fluid in pres- 
ence of a pycnocline, showing the robustness of dead-water 
phenomenon. We observe large amplitude nonlinear internal 
waves which are coupled to the boat dynamics, and we em- 
phasize that the modeling of the wave-induced drag requires 
more analysis, taking into account nonlinear effects. 



Dedicated to Fridtj of Hansen born 150 years ago (10 Oc- 
tober 1861). 



1 Introduction 

For sailors, the dead-water phenomenon is a well-known pe- 
culiar phenomenon, when a boat evolving on a two-layer 
fluid feels an extra drag due to waves being generated at the 
interface between the two layers whereas the free surface re- 
mains still. Interestingly, one finds reports of similar phe- 
nomena in the Latin literature when Tacitus described a flat 
sea on which one could not row a boat. North of Scotland 
and of Ge rmany, in the Agricola (ITacitusl l98ah and in the 
Germania (lTacitusll98bl) . 

This effect is only observed when the upper part of the 
fluid is composed of layers of different densities, due to vari- 
ations in salt concentration or temperature. An important 
loss of steering power and speed is experienced by the boat, 
which can even undergo an oscillatory motion when the mo- 
tors are stopped. 

In this paper, we present detailed experimental results 
on the dead-water phenomenon as shown in the video 



by IVasseur et al.l (l2008|) . The material is organized as fol- 
lows. In the remaining of this section, we briefly review the 
different studies of this phenomenon, either directly related 
to Ekman's work or only partially connected to it. Section |2] 
presents the experimental set-up. The case of a two-layer 
fluid is addressed in Sec.[3l followed by the case with a three- 
layer fluid in Sec. HI The more realistic stratification with a 
pycnocline above a linearly stratified fluid is finally discussed 
in Sec. [3 Our conclusions, and suggestions for future work 
are presented in section O 

1.1 Ekman's PhD Thesis 

V. W. Ekman was the first researcher to study in de- 
tail the origin of th e dead-water p henomenon. His work 
as a PhD student ( Ekmanl Il904 ) was motivated by the 
well-documented report from the Norvegian Artie explorer 
Fridtjof Nansen who experienced it while sailing on the Fram 
near "Nordeskiold" islands in 1893 dNanseni [18971) . 

Several aspects of the phenomenon have been described 
by Ekman, who did experiments with different types of boat 
evolving on a two-layer fluid. We note pi and hi the density 
and depth of the upper layer, and p2 and h2 those of the lower 
layer. 

i) First of all, the drag experienced by the boat evolving on 
the stratified fluid is much stronger than the one associated 
with an homogeneous fluid. This difference is due to wave 
generation at the interface between the two layers of fluid, 
pumping energy from the boat. This effect is the strongest 
when the boat's s peed is smaller than the maximum wave 
speed (lGiiill982l) . defined as 



p2-pi hih2 



(1) 
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as can be seen in Fig. [T] where the value of is indicated 
by the circle. A typical evolution of the drag versus speed 
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Fig. 1. Drag- speed relations given by Ekman for a boat dragged on 
a layer of fresh water (pi = 1.000 g cm^) of depth hi — 2.0 cm, 
resting above a salt layer {p2 = 1.030 g cm^), and compared to the 
homogeneous case "deep water" with a water level of ® 23 cm, 
and "shallow water" corresponding t o a smaller dep th d) 5 cm, and 
(D 2.5 cm (taken from (Fig. 8, PI. VI lEkmanL [19041) ). Experimental 
points are crosses and continuous lines are models. 



obtained by Ekman is shown in this figure. The experimen- 
tal points are crosses and are compared to linear theory of 
viscous drag in steady state (continuous lines). The local 
maximum of the drag for the stratified case is reached for a 
speed slightly smaller than c^. At high speeds compared to 
c^, the drag is similar to the quadratic law for homogeneous 
fluids, dominated by viscous drag. 

This behavior has bee n reproduced since t he n, and similar 



results were obtained by iMiloh et al.l (Il993h or IVosper et al, 
(Il999) for instance. 

ii) Another contribution of Ekman concerns the descrip- 
tion of the interfacial waves generated at the rear of the boat. 
Two types of waves, transverse and divergent, could be ob- 
served. It must be noted that only interfacial waves can be 
easily observed since the surface waves associated are of very 
weak amplitude, their amplitude being related to the interfa- 
cial waves ones with a ratio equals to (pa — pi) / {p2 — pi) — 
1/500 (where pa is the density of air). The free surface re- 
mains still at the laboratory scale. 

These waves are generated by a depression that develops 
at the rear of the boat while moving. The transverse waves 
can reach large amplitudes up to breaking. But they tend 
to disappear when the speed of the boat is greater than c^, 
where only divergent waves remain. 

Visualizations of these waves from numerical sim- 
ulations can 



also be 



waves 
observed in 



( Miloh etal.L 119931 



Yeung and Nguyeni Il999|) . Their properties are set by the 
Froude number only, which compares the mean speed of the 
moving object to c^. This will help us in our description 
later on. 

iii) Furthermore, a solitary wave at the bow of the boat can 
also be observed. This structure, a spatially localized bump 
is reminiscent of solutions of the Korteweg-DeVries (KdV) 
equation: it can evolve freely, conserving its shape when the 
boat is stopped. Otherwise, it remains trapped to the boat. 



iv) From these observations, Ekman gave an interpretation 
of the hysteretic behavior the speed of a boat can experience 
when evolving in dead-water. The analysis can be explained 
with the help of Fig. [T] showing the drag- speed relation es- 
tablished by Ekman from linear theory. If the force moving 
the boat increases such that the boat accelerates from rest, 
the speed of the boat can jump from 6 to 15 cm.s~^ when 
the boat overcomes the maximum drag. Similarly, when the 
force diminishes such that the boat decelerates from a speed 
larger than c^, a sudden decrease from 11 to 4 cm.s~^ oc- 
curs. The range of values between 6 and 11 cm.s~^ is thus an 
unstable branch inaccessible to the system, where dead-water 
phenomenon occurs. 

However it is important to emphasize that this analysis 
implies changing the moving force, hence not imposing a 
constant one. This observation is different from another re- 
mark made by Ekman which relates explicitly to the appar- 
ent un steady behavio r of the boat while towed by a constant 
force (lEkman[ri904 p. 67). He observed oscillations of the 
speed of the boat that could be of large amplitude compared 
to the mean value. He noticed that they occur when the boat 
is evolving at speeds smaller than c^, while amplitudes and 
periods of the oscillations depend on the towing force and the 
properties of the stratification. We emphasize here that this 
last property is not included in analytical approaches con 



Miloh et al 



sidering linear waves, such as ( Ekmani 1 190' 

1993l:lYeung and Nguyeni [l999 ). and seems to be an impor- 
tant characteristic of the dead-water phenomenon. 



1.2 Other dead-water related works 



Hughes and GrantI (Il978h took advantage of the dead-water 
effect to study the effects of interfacial waves on wind in- 
duced surface waves. The study relates the statistical prop- 
erties of surface waves to the currents induced by internal 

w aves. 

Maas and van HarenI (l2006h investigated if the dead-water 



effect could also be experienced by swimmers in a thermally 
stratified pool, offering a plausible explanation for unex- 
plained drownings of experienced swimmers in lakes during 
the summer season, but found no effect. One can argue that 
the stratification considered might have been inadequate for 
swimmers to generate waves and led to mixing of the ther- 
mocline mainly. An energ etic budget is given in a more de- 
tailed and idealized study ( Ganzevles et al.[ l2009|) where the 
authors also observed some retarding effects on the swim- 
mers. 



In a slightly different perspective, iNicolaou et al.l (Il995h 
demonstrated that an object accelerating in a stratified fluid 
generates oblique and transverse internal waves, the lat- 
ter can be decomposed as a sum of baroclinic modes with 
the lowest mode always present. Shishkina (2QQZ) further 
showed through experiments that in such a dynamical evo- 
lution, the baroclinic modes generated propagate indepen- 
dently of each other, although nonlinear effects must become 
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important when the ampHtude of the internal waves is in- 
creasing. 

1.3 Steady motion, body moving at constant speed 

Finally, numerous studies were focused on bodies evolving 
at constant speed within a stratified fluid. Some results can 
help the understanding of the dead-water experiments, espe- 
cially the internal waves at the rear of the boat and the wave- 
induced drag on the boat. 

In the case of a two-layer fluid, the drag on the boat 
is maximal when the Froude number, defined as the ra- 
tio of the boat speed U to the maximum wave speed 
given in Eg. ([T]), is slightly less than 1 ( Miloh et aP. 



19931: iMotygin and Kuznetsovl. 119970 : this is the subcrit- 
ical regime. The structure of the internal waves gen- 
erated and their coupling with surface waves confirm 
that the dead-water regime i s due to baroclinic waves 
only dYeung and Nguyen[ll999 ). 



These results obtained in a linear case c an be extende d 
when considering weakly nonlinear effects ( Bainesl Il995k . 
Fully nonlinear calculations are needed when the amplitude 
of th e waves reaches about 0. 4 times the depth of the thinner 
laver lGrue et"aD d 19971 1 1999b . 

In the case of a linearly stratified fluid, the drag 
is again maxima l for slight l y sub critical values of the 
Froude number fcree nsladel 12000 ). Nevertheless, the 
internal waves emitted can be very different depend- 
ing on the regime con sidered for the Froude num- 
ber ( Chomaz et al.[ Il993|). the loca ti on of the object be- 
ing at the surface dRottman et al.l 20041) or full y im- 



mersed (IHopfinger et al 



[ml 



IMeunier et al.. 120061). It is 



interesting to emphasize that in experiments done at constant 
speed, difficulties are often encounter ed to reach a steady 
state, as noticed by lVosper et al.l d 1999k for instance. 

Nevertheless, the dead-water phenomenon does not corre- 
spond to a constant speed evolution. By imposing a constant 
force to move the boat, or using a motor at constant power, 
the speed of the boat is free to evolve. The dynamical study 
of the problem is thus much richer than in the steady state 
case. 



2 Experimental setup 

Nowadays technologies give us the opportunity to gain more 
insight into the interactions between the interfacial waves and 
the boat, and improve our knowledge of the dead-water phe- 
nomenon. The experimental setup is described in Figure |2] 

We drag a plastic Playmobil© boat of width 10 cm (and 
with a fisherman and fishes to modulate its weight) with a 
falling weight in a 3-m long plexiglass tank of width 10.5 cm 
filled with a stratified fluid. A belt with fixed tension is used 
to move the boat with a constant horizontal force. The ten- 
sion of the belt is set to a constant value throughout all the 



experiments. The falling weight of mass m fixed to the belt 
sets the boat into motion, such that the constant force used is 
gravity mg. One must notice that the weights used are paper 
clips of few milligrams, since a very small force is required to 
propel the boat within the interesting regime. A magnet fixed 
to the boat is used to release it in a systematic way thanks to 
an electro-magnet outside the tank. 




Fig. 2. Experimental setup of a boat dragged by a falling weight. 
The vertical force is converted horizontally through pulleys and an 
horizontal belt of fixed tension. 



2.1 Stratification 

Different types of stratification have been used. The density 
profiles presented here are obtained using a conductivity and 
temperature probe (CT-probe) from PME ©. 

1) The two-layer fluid is composed of a layer of fresh wa- 
ter (density pi) colored with red food dye resting above a 
transparent layer of salt water (density p2> pi)- By siphon- 
ing the interface, the density jump extension can be reduced 
to a few millimeters but through successive experiments, dif- 
fusion and mixing make it widen with time, up to a few cen- 
timeters as can be seen in Fig.O 



^ -0.08 




/>(kg.m ^) 



N (rad.s-^) 



Fig. 3. Density profile p and Brunt- Vaisala frequency N for a two- 
layer fluid after a series of experiments. The density jump is ap- 
proximately 4 cm wide as observed in the zoomed window close to 
the free surface. 

ii) The three-layer fluid is obtained similarly, by adding an- 
other layer of salt fluid (density P3> P2) colored with green 
food dye from the bottom. 

iii) Finally, a continuously stratified fluid composed of an 
homogeneous layer above a linearly stratified part can be ob- 
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tained from adding a fresh laye r after fillin g the tank with the 
classic "two-bucket" method dHill l2002h . An example of 
such stratification with pycnocline is shown in Fig. (H along 
with the Brunt- Vaisala frequency associated with it. 



(a) 
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Fig. 4. Density profile p and Brunt- Vaisala frequency N obtained 
experimentally in the case of a linear stratification with a pycno- 
cline. 



2.2 Techniques of image analysis 

We record the dynamics of the system (waves -\- boat) using 
a black and white camera. Depending on the stratification 
considered, two different techniques are used to extract the 
dynamics. For the two and three-layer cases, the different 
layers are identified using food dye and the position of the in- 
terfaces r]{x^t) are extracted from highly contrasted images, 
as can be seen in Fig. [Jl The boat position with time x{t) 
and the free surface evolution are also obtained from these 
images. Several hypothesis must be highlighted when deal- 
ing with information obtained with this technique, that we 
called technique 1 : 

- we consider the interface between two layers as in- 
finitely thin, 

- we neglect diffusion of salt and dye, which is a slow 
phenomenon compared to the experiments (typically a 
day versus a few minutes), 

- we neglect the small scale evolutions of the interface, 
more specifically mixing events. 

The second te chnique used , technique 2 is based on syn- 
thetic schlieren (iDalziel et al. . 120001 120071) and can be real- 
ized simultaneously with technique 1. By computing corre- 
lations between images in the stratified case and a reference 
image with homogeneous water, one can measure the com- 
plete density gradient and not only the one associated with 
the internal waves field. In the domain considered, one can 
thus SiCCQSS dx p{x J z ^t) Sind dzp{x^z^t). In the absence of any 
motion, the second quantities have been integrated to mea- 
sure the density profile and it was compared successfully to 
the one obtained with the CT-probe. Allowing to quantify the 
large oscillations of strong density gradient, this technique 





Fig. 5. (a) Gray scale image during an experiment, converted into a 
(b) two-level black and white image with technique 1 . 



has been especially considered in the case of a continuous 
stratification with a pycnocline. 

We have checked that technique 1 gives a good indication 
of the position of the interface by recording synchronously a 
two-layer case with both techniques. As can be seen in Fig. [6] 
where the image (a) is recorded with technique 1 and im- 
age (b) with technique 2, the position of the interface (black 
line) extracted from the image (a) follows the main evolution 
of the density jump in (b) where the vertical density gradient 
is the strongest. 



& -0. 
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-0.05 0.05 



0.25 0.3 0.35 
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Fig. 6. Technique 2 for a two-layer fluid experiment: (a) color im- 
age (taken with camera) and (b) gray scale image over which is 
superposed the vertical gradient of density dzping cm~^, obtained 
using synthetic schlieren with a reference image with homogeneous 
water. The black line corresponds to the interface position extracted 
from image (a) and added to image (b) while the pink dashed line is 
the free surface. 



Although technique 2 can be used for all experiments, we 
will use technique 1 to analyze the two and three-layer cases 
since its computational cost is much less. 
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3 Revisiting Ekman's work 

Our first experimental investigation is dedicated to the two- 
layer case, in order to reproduce the results obtained by Ek- 
man. The main parameters are summarized in Table [T] 

3.1 Dynamics of the boat 

As can be seen in Table [H we have used several sets of pa- 
rameters in order to verify that the dead-water phenomenon 
is indeed a function of the relative density jump Ap = 
{p2 — pi)/p, the fresh layer depth hi compared to the salted 
one /i2 along with the waterline on the boat (changing 
hi) implies a different immersed cross-section and thus a dif- 
ferent viscous drag). In the following, we will refer to the 
immersed cross-section Sh for varying geometry of the boat, 
and which corresponds to w^xhi). 

For instance, the dead-water phenomenon is easy to repro- 
duce with the following parameters: hi=b cm, /12 = 14 cm, 
hb 2.4 cm, Ap = 0.0244 and the falling weight corre- 
sponding to a force Ft = 16.3 mN (also studied in Sec. 13.21 
Fig- [IS- The time evolution of the boat x{t) and its corre- 
sponding oscillating speed v{t) =dx/dt normalized by 
associated with this experiment, are shown in blue in Fig.H 
A similar case but with a force = 18.8 mN (also studied 
in Sec. 13.21 Fig.fTTI) is shown in green in Fig. [71 




t(s) 



Fig. 7. (a) Position of the boat x(t) (in m) and (b) its speed com- 
pared to the maximum phase speed v{t)/c^ for experiments re- 
vealing the dead-water phenomenon (blue and green) and for an 
experiment without oscillation (red). 

Ekman already commented on the fact that the speed can 
have large fluctuations \vmax — ^mm| ^ 0.027 m.s~^ com- 
pared to its mean value <v>~0.031 m s~^ (blue case), 
referring to the capricious nature of the phenomenon. Oscil- 
lations of the order of 85% of the mean value make it difficult 
to interpret the dead-water regime as stationary. The force 
imposed to move the boat being constant, such a temporal 
evolution implies a time-evolving drag due to the interfacial 
waves, and that is due to generation and breaking of waves 









* 








* > 



































0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 
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Fig. 8. Speed-Force diagram for three stratifications with a fresh 
layer h\ being respectively (*)0.06 m, (o)0.05 m and (<)0.04 m; 
with a salt layer /i2 = 0.14 m and Ap = 0.0247, for a boat con- 
figuration being 5'b = 12 cm^. For the unsteady regime, horizontal 
lines are drawn corresponding to the range of values of the speed 
in which it fluctuates. The continuous line corresponds to a homo- 
geneous case of depth 0.15 m. The vertical dashed lines gives the 
values of associated with h\, (dotted line) 0.10 m s~^, (dash- 
dotted line) 0.095 m s"^ and (dashed line) 0.087 m s"\ 

at the rear of the boat. This will be confirmed later on when 
we present the spatio-temporal diagram associated with the 
interface position. 

We also show in red in Fig.|7]the time evolution of x(i) and 
v{t)/c^ for the same parameters except Ft = 21 mN (also 
studied in Sec. 13.21 Fig. [TT]). Although only reached at the 
end of the recording, the boat tends to a steady evolution at 
constant speed, V£~0.14ms~^, after a transient state (for 
t < 15 s). The comparison between the red and green case 
also reveals the sensitivity of the system regarding the drag 
force. The temporal evolutions are initially very similar (for 
the first five seconds) before behaving differently. For suf- 
ficiently large speeds, the hypothesis of constant velocity of 
the boat seems reliable and can lead to a good interpretation 
of the evolution of the boat. 

Analogously to Ekman's results presented in Fig. [T] and 
in order to compare our results with his, we summarize in 
Figs. [8] and. [9] several time evolutions for the different strati- 
fications and boat configurations considered in a speed-force 
diagram. However, two comments must be made with re- 
spect to Ekman's presentation of his experimental results. 
When the boat is in an oscillating regime, Ekman takes the 
mean value as its speed and considers that the drag equals 
the moving force. This simplified representation of the data 
leads to cleaner diagrams although some strong hypotheses 
are hidden. 

In our case, we consider the moving force Ft instead of 
the drag force since both forces are equal only when the boat 
evolution is steady. We represent the speed by its mean value 
(symbols) and the range of values (lines) in which it fluctu- 
ates. 

Figure [8] combines the speed-force diagram associated 
with three different stratifications and one boat configuration 
{Sh = 12 cm^), whereas Fig. [9] gathers three boat configura- 
tions and one stratification. 
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Parameters 




symbols 


values 


units 


Tank 


dimensions 


LxHxW 


300 X 50 X 10.5 


cm^ 


Traction 


belt tension 


T 


-0.35 


N 


force 


Ft 


0.011-0.035 


N 




dimensions 


Tji X X iiiu 


20.0 X 5.0 X 10.0 


cm*^ 


Boat 


immersed section 




12.0-24.0 


cm^ 




mass 


M 


171-343 


g 


Fluid 1 


density 
depth 


Pi 
hi 


0.999-1.005 
2.0-5.0 


gem 
cm 


Fluid 2 


density 
depth 


P2 

h2 


1.010-1.030 
5.0-15.0 


gem 
cm 




mean density 




-1.01 


gem 




density jump 


Ap- P^-P^ 


0.01-0.1 






maximum phase speed 




0.03-0.2 


m s~ 




Froude number 




0.2-2 






Reynolds number 




400 - 10000 





Table 1. Experimental parameters used for the experiments with a two -layer fluid. 




0.02 0.04 0.06 0.08 0.1 0.12 0.14 0.16 0.18 0.2 

vt (m s~-^) 

Fig. 9. Speed-Force diagram for three boat configurations, (*) — 
Sh = 12, (<) -Sb = 18, and (o) - 5'b = 24 cm^ in a stratification 
with hi = 0.06 m, /i2 = 0.10 m and Ap = 0.0145. For the unsteady 
regime, horizontal lines are drawn corresponding to the range of 
values of the speed in which it fluctuates. The continuous line cor- 
responds to a homogeneous case of depth 0.15 m with Sb = lS cm^. 
The vertical dashed line gives the value of 



change of the moving force leads to a completely different 
evolution of the boat. One can thus argue there exists a range 
of "unattainable" values for the velocity of the boat, as dis- 
cussed in Sec. 11.11 § iv. However, since the moving force 
is kept constant in our case, this aspect is different from Ek- 
man's presentation of the phenomenon. Similar comparisons 
have been made for the influence of the density jump Ap and 
the salt layer depth h2 but are not shown here for the sake of 
concision. 

Finally, several general comments are in order 

- the wave-induced drag is always stronger than the one 
in homogeneous water (viscous drag only), 

- a regime of low velocities exists with a non constant 
drag and large fluctuations in the boat speed (no steady 
state can be reached). 



These two diagrams are in agreement with the established 
results by Ekman that the wave-induced drag increases along 
with the fresh water depth hi when a steady state exists 
(Fig. [5]), and depends also on the boat geometry (Fig. O. 
As the immersed section of the boat is larger, the drag gets 
stronger. The best coupling occurs for hb — hi/2 where the 
wave induced drag is maximum and the boat is not deep 
enough to generate mixing at the interface. Furthermore, 
both diagrams show that when the speed of the boat is os- 
cillating, it never reaches values larger than c^. It is impor- 
tant to notice that the dynamical regime where the velocity 
of the boat is closed to cT is difficult to scan since the tiniest 



- this regime is always associated with velocities smaller 
thanc^. 

From these remarks, we can discriminate the two regimes de- 
scribed previously according to the Froude number Fr which 
is defined as the ratio of the speed of the boat (that can be 
its mean value <v> or its limit V£) to the maximum phase 
speed of the interfacial waves. It is important to realize here 
that the Froude number is intrinsically a time evolving pa- 
rameter, but we will use characteristic values to identify the 
regime considered. 

We now focus on the waves according to this classifica- 
tion. 



?7(x,t) (m) 




X (m) 



Fig. 10. Ballistic regime (Fr > 1). Spatio-temporal diagram of 
interfacial displacements r]{x,t) (in m). The solid black lines are the 
bow and stem of the boat. Experimental parameters: hi = 5.0 cm, 
h2 = 14.0 cm, pi = 0.9980 g cm"^ p2 = 1.0227 g cm"^ Sb = 
24 cm^F^ = 21mN. 
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?7(x,t) (m) 
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Fig. 11. Oscillating regime (Fr < 1). Spatio-temporal diagram of 
interfacial displacements r]{x,t) (in m). The solid black lines are the 
bow and stem of the boat. Experimental parameters: hi = 5.0 cm, 
h2 14.0 cm, pi 0.9980 g cm"^ p2 = 1.0227 g cm"^ Sb 
24 cm^ Ft = 18.8 mN. 



3.2 Interfacial waves 

We present here observations of the interfacial position 
r]{x^t) generated by the boat in two different frames, the lab- 
oratory reference frame and the one moving with the boat. 
The latter is not Galilean since the boat does not evolve at 
constant speed, but offers a better understanding of the cou- 
pled dynamics of the waves with the boat. 

Visualizations of the interfacial positions are represented 
in spatio-temporal diagrams where elevations of the interface 
are in red and depressions in blue (see Figs. [TOl and ITT] for 
instance). The bow and stem of the boat are superimposed 
on these diagrams (as black lines). 

Depending on the Froude number being larger or smaller 
than 1, the wave dynamics is very different. We illustrate 
each case with an example corresponding to a stratification 
with hi = 5.0 cm, h2 = 14.0 cm, pi = 0.9980 g cm"^, p2 
1.0227 g cm-^, leading to = 0.094 m s-^ The boat used 
corresponds to an immersed section = 24 cm^ . 

3.2.1 In the laboratory frame 

Let us first consider the case where a steady state is reached, 
illustrated in Fig.[TOl With a moving force Ft = 21 mN, the 
limit speed reached by the boat is V£ 0.15 m s~^ which 
gives Fr = 1.6. In this supercritical case, one can observe 
that the boat escapes from the wave train it generates at start- 
up. A depression (the blue region) remains attached below 
the boat and propagates with it at the speed of the boat. For 
times larger than 40 s, the boat meets with the end of the tank 
and stops (around x = 2.5 m). The depression below it is thus 
expelled and propagates in the other direction after reflecting 
on the side of the tank. This case is in the ballistic regime. 



The oscillating regime, more atypical, with oscillations of 
the boat is illustrated in Fig. [TT] and corresponds to Ft = 
18.8 mN. The mean speed of the boat is < v 0.035 ms~^, 
which gives a Froude number of 0.4. As the boat evolves, the 
amplitude of the interfacial waves grows up to a maximum 
value where the waves catch up with the boat and break on 
its hull, around t = 20 s. This corresponds to the moment 
where the drag force is maximal since the boat is almost 
standing still. The vessel can speed up again after the first 
wave breaks and the same dynamics is reproduced. As noted 
earlier, a depression located below the boat evolves with the 
boat when it is accelerating. It evolves freely after the boat is 
stopped (t > 20 s) and propagates to the right first and to the 
left for t > 50 s after reflection at the side of the tank. This 
depression even bounces a second time on the left side of the 
tank and its shape remains almost unchanged (t > 75 s). Al- 
though not presented here, we have verified that this solitary 
wave propagating upstream is well described by a Korteweg- 
de Vrie s model since its amplitude is always smaller than 



0.4 hi (iGrue et al.L 119991) . When the boat starts again, the 



process is repeated and a new depression is generated. 
3.2.2 In the frame of the boat 

More information can be obtained when following the dy- 
namics in the frame associated with the boat. We are more 
specifically interested in what sets the amplitude and fre- 
quency of the oscillations of the boat when Fr <1. 

In Fig. [T2I we have superposed the interfacial displace- 
ments in the frame of the boat (represented as a black rect- 
angle) for different times extracted from Fig. [TT] We can 
observe that the waves get closer to the stem of the boat as 
their amplitude grows. This is mainly due to the decrease of 
speed of the boat since we can see in Fig. [TT] that the wave 
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Xh (m) 

Fig. 12. Fr <1. Interface evolution rj{x,t) (in m) from its rest 
position hi for several times during the slow down of the boat, in 
the frame of the boat which is represented by the black rectangular 
box. Same parameters as in Fig.fTTI 



crests evolve with a constant speed. The nonlinear nature of 
the interfacial waves is also clearly visible as the wave fronts 
steepen, the limit of growth being set here by the time when 
the waves hit the hull (27 < t < 32 s). These downstream fea- 
tures are fully nonlinear (amplitudes larger than 0.6 hi) and 
have similar a spects to observations made in the transcriti- 
cal regime by (Mel ville and Helfrichl [TqStI) and reproduced 
numerically by (G rue et al.l Il997k . although corresponding 
here to waves in a non-galilean frame. It is moreover sur- 
prising since as observed on the green curve in Fig. |7] the 
Froude number does not exceed 0.7. These observations add 
to the comments made previously concerning the inadequacy 
of models based on steady object in stratified fluid. Here 
again we see upstream the solitary depression below the boat 
escaping at the bow as a single oscillation. 

For times t > 32 s, one can observe that the decapitated 
first crest disappears below the boat, the next wave depres- 
sion takes position below the boat and the initially second el- 
evation of the packet becomes the closest to the stern. Hence, 
during an oscillation of the boat, the wave train slides of one 
wavelength compared to the boat, repeating itself. No change 
in the dominant wavelength of the wave train is visible. 

We can now estimate the frequency of the oscillations of 
the boat to be 



h- 



X 



(2) 



with Cg the group velocity of the wave train at the rear of the 
boat and A its dominant wavelength, corresponding approxi- 
mately to the distance between the first two wave crests. 

The experiment presented in Fig. [TT] allows an observation 
of the steepening process but does not present many oscilla- 
tions. In order to confirm the validity of formula (|2), we con- 
sider another experiment with = 16.3 mN in Fig. [13] From 
Fig.[T3](a), one can estimate the frequency of the oscillations 
of the boat to be ~ 0.07 — 0.08 Hz. The group speed and 
wavelength of the wave train obtained from the diagram in 
Fig. [13] (b), corresponds to Cg - 0.035 - 0.040 m s"^ and 
A ~ 0.5 m. Indeed, we obtain Cg/X^ 
ing very good agreement with Z^. 




Fig. 13. Oscillating regime {Fr < 1). Panel (a) presents the speed 
of the boat versus time, while panels (b) and (c) shows the spatio- 
temporal diagrams of the interfacial displacements in the laboratory 
frame and in the one of the boat respectively. The black lines rep- 
resent the bow and stem. Experimental parameters: /ii = 5.0 cm, 



/i2 14.0 cm, pi 0.9980 g cm" 
24 cm^ Ft = 16.3 mN. 



P2 = 1.0157 g cm-^ Sb-- 



0.07-0.08 Hz show- 



After revisiting experiments similar to the ones performed 
by Ekman, we have considered new situations. Indeed, dead- 
water phenomenon can be studied in more complex (and real- 
istic) stratifications, when the fluid has more than two layers 
of different densities. We now turn to new experiments with 
a three-layer fluid or with a continuously stratified fluid with 
a pycnocline, where the complex dynamics has also been ob- 
served. 



4 Experiments with a three-layer fluid 

We consider a three-layer fluid of depth hi and density pi, 
i = 1,2,3. Two interfaces must be considered now, r]ij{x,t) 
corresponding to the interfacial displacements between lay- 
ers i and j. A picture summarizes the setup in Fig. [14] and 
the main parameters are given in Table |2] Only one given 
stratification will be discussed here. 

We use technique 1 again to extract the positions of the 
two interfaces r]i2 et 7/23, along with the position of the boat. 

4.1 Analytics 

Although the three-layer case is a simple extension of the 
well documented two-layer one discussed in previous sec- 
tion, we here reproduce explicit expressions for the phase 
speed of the two eigenmode solutions of this problem, 
w hich can also be find es tab lished in the most general forn i 
bv lRusas and Grud (l2002h or lXiao-Gang and Jin-Baol (l2006h 
for example. 
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Parameters 



symbols 



values 



umts 



Fluid 1 



density 
depth 



hi 



0.9967 
5.0 



gem 
cm 



Fluid 2 



Fluid 3 



Interface 7712 



Interface 7723 



Modes s/a 



density 
depth 

density 
depth 



P2 
P3 



2 

P2-P1 
P12 



mean density P12 = 

density jump Api2 

maximum phase speed c^,i2 = \J ^pi2g 

Froude number Fr 12 — 

mean density p23 - 

density jump Ap23 

maximum phase speed 

Froude number Fr2 

maximum phase speed 

Froude number Fr, 



1^2 



hi+h2 



f ,12 
P3+P2 
2 

P23 



/,,23-yAp23^7;;^ip7^ 



s/a " 



^(/),23 
s/a 



1.0079 
3.0 

1.0201 
5.5 

1.0023 
0.0112 

0.0147 
4.1-7.5 

1.014 
0.012 

0.0150 
4.0-7.3 
0.073/0.044 
1.5/2.4-0.6/0.9 



Table 2. Experimental parameters used for the experiments with a three-layer fluid. 



gem 
cm 

gem- 
cm 

gem" 



gem 




Fig. 14. Picture illustrating the dead-water phenomenon for a three -layer fluid. 



Pi 



tanh(/c/ii) tanh(/c/i^+i) 



Perturbations of a three-layer fluid can be described by two 
types of harmonic interfacial waves, verifying the dispersion 
relation expressed as the following determinant 

with, for i = 1 , 2 

c^i = 9\Pi- Pi+i)^- 

/3 = . 

/csinh(/c/i2) 

The fourth order polynomial in associated with ([3]) leads 
to two types of interfacial waves. Using the notations in- 
troduced in Table O and considering the long-wave limit 
(tanh(/c/ii) ^ sYnh{khi) khi) along with the limit of small 
density jumps (Apij <C Pi^Pj and pi ~ p^ij), we can ex- 
press the two phase velocities associated with each solution 

~ 1 1/2 



(4) 



(5) 



^1 + ^2 + ^3 



,(6) 



Cn = 



with 



^ AP1 2/11 (/i2 + hs) + Ap23h3{hi + /12) - A 
^1 + ^2 + ^3 



1/2 



,(7) 



A = ^ [Ap?2^?(/i2 -^hsf^ Apl^hlihi + /i2)^ 



-2Api2Ap23hih2hs{hi^h2^hs - 



hihs 
h2 



1/2 



. (8) 



These two phase speeds are associated with symmetric and 
anti- symmetric oscillations of the interfaces, 

Vs{x,t) = i(7/i2(x,t)+7?23(^,^)), (9) 

Va{x,t) = ]^(r]i2{x,t)-r]2z{x,t)). (10) 

We will refer to mode- 5 and mode-a for r]s and r]a respec- 
tively. We also define two Froude numbers Fr^ia = ^/c^a 
which are the appropriate non-dimensional numbers to de- 
limitate the oscillating regime as we will see in the follow- 
ing. Numerical values from the experiments are summarized 
in Table O 
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r]s{x,t) (m) 
0.03 



0.02 




0.01 



0.02 



0.03 



?7a(x,t) (m) 
0.015 



0.01 



0.005 




■0.005 



0.01 



■□.□15 



X (m) 

Fig. 15. Ballistic regime (Fvs > 1). Spatio-temporal diagram of interfacial displacements r]s{x,t) and r]a{x,t) (in m). The solid black lines 
are the bow and stem of the boat. The magenta (resp. white) dotted line is an indication of the speed of propagation of 0.028 m s~^ (resp. 
0.0575 m s~^). Experimental parameters are given in Table|2]and Sb — 24 cm^, Ft — 23.5 mN. 



To emphasize that the two interfaces cannot be seen as in- 
dependent oscillatory structures, since they are coupled by 
the intermediate layer, we give in table |2] the values of the 
Froude numbers for the two-layer equivalent solutions, Fr 12 
and Fr23. The values associated with the experiments do not 
allow the sub/supercritical classification. 

4.2 Case Fvs > 1 and Fra > 1 

We first consider the case where the vessel is going faster 
than the fastest wave. The experimental parameters associ- 
ated with the results shown in Fig. [15] are Sb = 24 cm^ and 
Ft = 23.5 mN. The limit speed of the boat is 0.11 m s~^, 
which gives Fts = 1.5 and Fra = 2.4. We only present the 
interfacial displacement as r]s and r]a because the spatio- 
temporal diagrams for 7^12 and 7^23 are very similar and do 
not clearly exhibit the link between the waves and the boat 
dynamics. 

Both modes are generated initially but the amplitude of 
mode-s is more than two times larger than that of mode- a 
which is at the limit of detection. Symmetric oscillations 
form a clear wave train whose front propagates at a speed of 
0.0575 m s~^ , while showing a similar spreading of the wave 
crests as in the two-layer case: this is an indication of disper- 
sion effects. Speed of propagation of mode- a is0.028ms~^. 

4.3 Case Fts < 1 and Fva ^ 1 

We now consider a subcritical case. It is however difficult to 
define clearly this domain with respect to both modes. The 
experimental parameters associated with the results shown in 
Fig.[T6lare Sb = 24 cm^ and Ft = 20.6 mN. The mean speed 



of the boat is <v >c^ 0.06 m s~^ which gives Frg = 0.6 and 
Fra = 0.9, with oscillations of the order of 0.04 m s~^. 

As observed previously, both modes are generated when 
the boat starts. Mode-s remains two times larger than mode- 
a but they are of noticeable amplitude. The symmetric mode 
evolves with the boat and reproduces the amplification and 
steepening processes observed in the two-layer case. The 
mode- 5 breaks on the boat when its amplitude is the largest, 
and a depression (symmetric for both interfaces) is expelled 
at the bow. The characteristics of the mode-s, wavelength 
Xs ^ 0.4 m and group velocity Cg^s — 0-06 m s~^ lead to an 
oscillating frequency of 0.15 Hz, too slow for this visualiza- 
tion. 

Concerning mode-a, its spatial structure is close to a soli- 
tary wave train, with a group velocity Cg^a — 0.0325 m s~^, 
close to (being 0.044 m s~^). Each time it is generated, a 
strong acceleration of the boat occurs. 

4.4 Discussion 

The unsteady behavior associated with dead-water is still ob- 
served in the three-layer fluid, with strong analogy to the 
two-layer case. The stratification considered being more 
complex, we must consider two baroclinic modes associated 
with symmetric and anti- symmetric oscillations of the inter- 
faces and that are also referred to as mode-1 and mode-2 in 
the literature. 



Rusas and Grud (l2002h present solutions of the nonlin- 



ear equations that have strong similarities with our obser- 
vations. More specifically, the spatio-temporal diagram of 
mode-a exhibit solitary waves of mode-2 with oscillatory 
short mode-1 waves superimposed. This is in agreement with 
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r]s{x,t) (m) 
0.03 



0.02 




0.01 



0.01 



0.02 



-0.03 



r]a{x,t) (m) 
0.015 



0.01 




0.005 



■0.005 



■0.01 



■□.□15 



Fig. 16. Oscillating regime (Fvs < 1). Spatio-temporal diagram of interfacial displacements r]s{x,t) and rja{x^t) (in m). The solid black 
lines are the bow and stem of the boat. The magenta (resp. white) dotted line is an indication of the speed of propagation of 0.0325 m s~^ 
(resp. 0.06 m s~^). Experimental parameters are given in in Table|2]and Sh — 24 cm^, Ft — 20.6 mN. 



close values of the experiments with the numerical calcula- 
tions (Boussinesq limit, hi/h^, — ! and hi/h2 — 1.7). 

Perturbations generated by the boat give birth to both 
modes, especially when the acceleration of the boat is impor- 



tant. T his result is consistent with the study ofiNicolaou et al,. 
(1995), verified experimentally by Ikobeyl d 1997b , stating that 



an accelerated object in a continuously stratified fluid, with 
a Brunt-Vaisala frequency N{z), excites a continuum of 
modes whose vertical profile w{z) described by 



dr'w 



N\z) 



l]w{z) = 0, 



(11) 



along with the boundary conditions for the vertical velocity 
to be zero at the top and bottom. 

Finally, we have observed that the mode-1 is strongly cou- 
pled to the dynamics of the boat, which corresponds to the 
fastest wave propagating in this stratification. A weak mode- 
2 is associated to noticeable acceleration of the boat, but 
evolves freely from the other mode. 



5 Continuously stratified fluid with a pycnocline 

We have observed that several waves are generated when the 
boat evolves in a complex stratification. In the case of a lin- 
early stratified fluid, an infinite number of modes can propa- 
gate. We have actually considered the case of a linearly strat- 
ified fluid with a pycnocline (see Fig. Several reasons 
can be invoked. It keeps the stratification undisrupted when 
the boat evolves in the top layer, it allows larger vertical dis- 
placements at the density jump leading to larger amplitudes. 



it can be modeled as a coupling between interfacial and in- 
ternal waves, and it corresponds to a more realistic setup in 
comparison with observations made in natural environment. 

Experimental parameters associated with the experiments 
presented are in Table [3] We use technique- 2 in order to ob- 
serve both interfacial oscillations and internal waves. 

Similarly to the three-layer case, the complex stratification 
considered here allows different definition for the Froude 
number. We know from the previous study that the one asso- 
ciated with the internal waves and based on the phase speed 
of the first vertical mode, Frh = Uir/Nh, is the adequate 
one although its value is actually not so different from the 
one obtained from the equivalent two-layer fluid. 

5.1 Modal decomposition 



As discussed in Sec. 14.41 we expect the non-stationary evolu- 
tion of the boat to generate several modes. We anticipate here 
by presenting the expected mode structures associated with 
the stratification (Fig.|4]) and obtained from solving Eq. (ITTI) . 
From the vertical structure of the vertical component of the 
velocity obtained from (fTTI) . we can deduce the correspond- 
ing density gradients, since 



UJ 

then 



(12) 



(13) 



where is a reference density uj the frequency of the mode 
and kx its horizontal wavenumber. 
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Parameters 






symbols 


values 


units 


Fluid 1 


density 
depth 




Pi 
hi 


0.9978 
5.0 


gem" 
cm 


Fluid 2 


density 
depth 

BV frequency 




N 


1.007-1.028 
20.0 
0.9 


gcm~ 
cm 

rad s~ 


Interface 7712 


mean density 
density jump 






1.0024 
0.0092 


gem" 


maximum phase 
Froude number 


speed 


Fri2 = ^ 


0.06 
0.6-2.3 


m s~^ 


Modes 


maximum phase 
Froude number 


speed 


Nh 

Fru — — 


0.072 
0.5-1.9 


m s~^ 



Table 3. Experimental parameters used for the experiments with a linearly stratified fluid with a pycnocline of depth h — hi + h2. 




dxpn 



dxpn 



Fig. 17. Normalized horizontal density gradient of the first three 
vertical modes of stratification in Fig. |4] with (a) cji = 0.25 rad s~^ 
and (b) uj2 = 1.1 rad s~^. 



Figure [T7] gives the vertical profiles (of the horizontal den- 
sity gradient) of the first three modes pn for two specific fre- 
quencies, (jOi = 0.25 rad s~^ in (a) and CJ2 = 1-1 rad s~^ in 
(b). The former corresponds to a mode propagating in the 
total depth of the fluid (since uji < N) whereas the latter is 
trapped in the pycnocline. (uj2 G [0.9; 1.8] rad s~^). 

The modal decomposition depends on the frequency of 
study. In our case, the dead-water phenomenon is a non- 
periodic state and we expect perturbations that are not har- 
monic modes. However, from previous observations, it is 
reasonable to expect nonlinear evolution of modes evolving 
freely from the others. Based on this assumption, we expect 
the following projection to catch the key features of the dy- 
namics. 

The computed modes in Fig. [T7] will help us identify dom- 
inant structures associated with specific frequency bands, but 
it is important not to apply too selective time filtering to the 
data in order to conserve information about their complex dy- 
namics. Thus we extract the mode amplitudes dn{x^t) after 
the two following steps done at each x-location: 

1. a large band filtering centered around uj with a band- 



width Auj, 

2. a projection of the vertical data onto the selected basis 
associated with the frequency band chosen uj 



1 

H 



H 



dxp'^{z).dxp{x,z,t)dz. 



(14) 



5.2 Case Frn > 1 



The first experiment considered corresponds to the experi- 
mental parameters in Table [3] with ^'5 = 12 cm^ and Ft = 
21.3 mN. The limit speed reached by the boat is Vi = 
0.14 m s~^ which corresponds to a supercritical case (Frh = 

I. 9). 

Figure [18] presents time series of the density gradient ex- 
tracted on a vertical cut at = 0.05 m where the origin is 
taken at the bow of the boat at t = 0. Interfacial displace- 
ments obtained from the evolution of the maximum of the 
vertical density gradient are given in white and are of less 
than 0.01 m. The fast evolution of the boat perturbs weakly 
the stratification and the signal associated with the internal 
waves field, for z G [-0.25,-0.05] m, is weak. 

Two types of internal waves can still be discriminated. Os- 
cillations of the full water depth are triggered when the boat 
passes by at t ^ 6 s. Later on, some turbulent fluctuations in 
the pycnocline lead to small scales forcing radiating waves 
downwards (t G [10; 20] s). The modal projection leads to no 
clear signal whatever the frequency considered and are not 
given here. 

5.3 Case Frh < 1 

The same configuration is used but with a force Ft = 

II. mN. The oscillatory dynamics of the boat around a 
mean speed <v>o:i 0.036 m s~^ corresponds to a subcrit- 
ical case with Frh = 0.5. 




"^^^0 5 10 15 20 25 30 95 40 45 50 

t(s) 

Fig. 18. Frh > 1. Time series of (a) dzp{xo,z,t) and (b) 

dxp{xo,z,t) extracted at xq = 0.5 m, in g/cm^. 




-0,25 ' ' ' 

Q 10 20 30 40 50 60 

t(s) 

Fig. 19. Frh < 1- Time series of (a) dzp{x,z,t) and (b) dxp{x,z,t) 
extracted at x = 0.67 m, in g/cm^. 



We represent the time series extracted at the location, 
X = 0.67 m (Fig. [191). The vertical cut corresponds to a re- 
gion close to where the boat will stop and its speed is al- 
ready decreasing. The stopping point of the boat is located at 
xc^l.bm and occurs at t 50 s. 

In Fig. \T9\ the displacements of the pycnocline is of the 
order of 0.01 m as in the supercritical case. Nevertheless, 
the internal waves field is much more intense and the ampli- 
tude of the waves generated are larger as the boat evolves. 
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(a) (b) 




X (m) X (m) 

Fig. 20. Frh < 1- Instantaneous horizontal density gradients 
dxp{x,z,to) at times (a) to = 32 s and (b) to = 61 s, in g/cm^. 



suggesting the same nonlinear evolution of the waves as for 
interfacial waves. 

The spatial structure seems to be dominated by a mode-1 
shape close to the boat, which is even more obvious in in- 
stantaneous visualizations of the density gradients as shown 
in Fig. [20] (a). Other modes are also present in the wave field 
but appear later in time and are hard to discriminate from the 
radiated waves generated by turbulent patches at the pycno- 
cline (Fig.[20](b)). 

In order to quantify in more detail the amplitude of each 
generated mode, we project the time series of the vertical 
structures extracted at different x-locations on the modal ba- 
sis described before and computed at five different frequen- 
cies, uj = [0.25,0.50,0.75,1.0,1.25] rad s-^ Large band fil- 
tering in time is done, with Auj = 0.25 rad s~^, so as to keep 
track of the nonlinear evolution of the amplitude of the modes 
with time. One can observe in Fig. [21] (a) that the ampli- 
tude of the mode-1 structure is indeed dominant since it is 
emitted for all frequencies considered, whereas the mode- 2 
in Fig. [21] (b) is of comparable amplitude at the lowest fre- 
quencies (cjO = 0.25 and 0.50 rad s~^) but is absent for larger 
values of uj. Higher modes (n > 3) are even weaker in am- 
plitudes and only present at low frequencies (Fig. [21] (c)). 

We notice that the initial perturbations for these low modes 
are generated below the boat since the black rectangles in 
Fig. [21] (a) and (b) are associated with a ramping amplitude 
of modes 1 and 2. The longer the rectangle is, the smaller 
the speed of the boat is, leading to larger initial perturbations 
which are very similar to what has been observed in the two- 
layer case (see Fig. [12]). 

Furthermore, the apparently cnoidal (sharp crests and flat 
troughs) oscillations associated with the mode-1 internal 
waves observed in Fig. [21] (a) as the modes propagate, are 
a signature of the nonlinear dynamics of the internal waves. 
Finally, by giving in all images the maximum phase speed of 
the waves with the tilted lines in Fig. [21] we can verify that 
the waves crests propagate at a smaller speed than 7.2 cm s~^ 
but it is difficult to extract a constant speed with propagation 
for each mode. 
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5.4 Discussion 




^=0.25 rad/s 

w =0.5 rad/s 

w=0.75 rad/s 

a;=l rad/s 

a;=1.25 rad/s 




10 20 30 40 50 60 70 

t(s) 

Fig. 21. Fr^ < 1. Amplitudes dn{xo,t) with time of the projected 
modal structure at several xo -locations for modes n=l (a), 2 (b) 
and 3 (c) of the horizontal density gradient in arbitrary units. Dif- 
ferent colors correspond to different frequencies associated with the 
modal basis. The tilted dashed lines corresponds to the charac- 
teristic speed of propagation c^^ = 7.2 cm/s. The black rectangle 
represents the boat passing through a vertical cross-section at each 
xo -location studied. 



The dead-water phenomenon has been observed in the more 
general case of a linearly stratified fluid with a pycnocline. 
The complex dynamics of the boat is coupled to the first 
mode of the stratification with frequencies below the max- 
imum value of N{z). 

Although the wave field associated with the subcritical 
regime is analogous to the previous layered stratifications 
considered, the supercritical regime is different in nature 
since it consists mainly of radiated waves from turbulent per- 
turbations in the pycnocline. This is similar to observations 
of the radiated wave field associated with an object movi ng 
at constant speed in a stratified fluid (IRottman et al .1120041) . 



6 Conclusions 

By revisiting the historical experiments of Ekman's PhD 
Thesis and extending it to more general stratifications, we 
have shown the robustness of the dead-water phenomenon. 
The experimental techniques used have revealed new insights 
on the century-old problem. 

One important characteristic of the dead-water phe- 
nomenon is that the dynamics of the boat is coupled to the 
fastest mode of the stratification considered, although sev- 
eral modes are generated at each acceleration of the boat. 
Furthermore, the nonlinear features of the phenomenon must 
be considered in order to describe analytically it s unstea dy 
nature. Cla ssical models such as presented by mad or 
Miloh et al.l are not sufficient. 

It would be of great interest to provide an analytical de- 
scription of the coupled dynamics of the boat and the waves 
it generates in order to take into account unsteady behaviors. 

An extended study of the problem in a three-dimensional 
setup is in preparation but we are still developing a collabo- 
ration with Playmobil ©. 
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